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Abstract 

^— ^ The phase-space path-integral approach to the damped harmonic oscillator is analyzed beyond the Markovian approximation. It 
^-^ is found that pairs of nonclassical trajectories contribute to the path-integral representation of the Wigner propagating function. 
p^J Due to the linearity of the problem, the sum coordinate of a pair still satisfies the classical equation of motion. Furthermore, it is 

K^, shown that the broadening of the Wigner propagating function of the damped oscillator arises due to the time-nonlocal interaction 

j^ mediated by the heat bath. 

f^ Keywords: harmonic oscillator, dissipation, non-Markovian, phase space 
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1. Introduction 
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In quantum chemistry, biophysics, and more recently in some 
problems in quantum computation, the time evolution of sys- 
tems of large size and at high excitations is in the center of 
interest. Generally, a full quantum treatment of such systems 
is practically impossible. In order to get a first insight into the 
O"^behavior of such systems, they are usually propagated using the 
'associated classical equations of motion. This approach, which 
^— I disregards all possible quantum effects, is known as molecular 
K>" 'dynamics (MD). Often, the scales involved in a problem call 
^^ -for a semiclassical approximation where molecular dynamics 
^yJ would constitute the zeroth order. To next order, one obtains 



o 
o 



>< 



*y^ the well-known van Vleck-Gutzwiller propagator ll], |2l] . This 
'quantity is arguably one of the most fundamental elements of 
semiclassical theory. It constitutes, e.g., the starting point for 
the derivation of the celebrated Gutzwiller trace formula yD and 
also for the semiclassical reaction rate theory ||4J,|5|]. 
■ Traditionally, semiclassical methods have been developed in 
"-H position representation. This implies, however, that i) determi- 
hants arising from the projection along momentum appear in 
^ the prefactor, ii) they are formulated in terms of double-sided 
■ " boundary condition problems and iii) are based on wave func- 
tions. Thus, they do not speak the same language as classical 
mechanics, which certainly would be desirable in order to de- 
scribe and study the quantum-classical transition i^]^- In view 
of these drawbacks, semiclassical methods in phase space offer 
a conceptually clearer and possibly a numerically more efficient 
approach i9l4llll. 

Recently, significant progress has been achieved as to semi- 
classical approximations to the Wigner propagator correspond- 
ing to the unitary quantum dynamical group. The case of 
damped nonunitary time evolution, giving rise to a semi-group, 
is different even in its basic mathematical structures and largely 
requires to reconsider the phase-space approach. In order to 
clearly distinguish the two cases, we shall use the term "prop- 



agating function" instead of "propagator" in the context of this 
semi-group. 

To be sure, a semiclassical phase-space analysis of open 
quantum systems has recently been presented in the Markovian 
limit [12 1 on the basis of the Lindblad master equation which 
can be extended to the general non-Markovian case il3ll for 
a system bilinearly coupled to a bath consisting of harmonic 
oscillators. In order to illustrate the latter approach by means 
of an analytically solvable case, we here present a study of 
the damped harmonic oscillator in phase space. While certain 
aspects will turn out to be specific for this linear system, we 
expect the overall picture to be useful to understand the more 
general case of nonUnear systems for which details will be dis- 
cussed elsewhere lll3ll . 

In an attempt to make the discussion accessible to readers 
with different backgrounds, we review in the next three sec- 
tions the basic ideas of the semiclassical phase-space approach, 
the influence functional theory, and the quantum damped har- 
monic oscillator as far as they are needed for the following. 
Section[5]will then contain the main results for the propagating 
function in phase space. We will argue that a Gaussian broaden- 
ing arises due to the nonlocal interaction of the system degree of 
freedom mediated by the coupling to the environment. Further- 
more, we demonstrate the nonclassical nature of the trajectories 
contributing within the phase-space path-integral representation 
of the Wigner propagating function. 

2. Unitary evolution of the density matrix 

Before discussing the case of a damped quantum system, it is 
useful to fix the main ideas and the notation by considering the 
unitary time evolution of an isolated system degree of freedom 
in real space as well as in phase space. In later sections, we will 
be forced to employ a density operator to describe the state of 
a damped quantum system. Therefore, we express the state of 
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our isolated system S, which could well be a pure state, also in 
terms of a density matrix denoted by ps in the following. 

For an isolated system described by a time-independent 
Hamiltonian H5, the time evolution of an initial density matrix 
ps(0) is determined through the unitary time evolution operator 



(1) 



f/(f) = exp(--//sr 



and its adjoint operator C/ ' (f) by means of the relation 

ps(r) = U{t)ps(OWHt). 
In position representation, this expression turns into 



(2) 



Psiq'-l, q'-, t) = I dq'^dq'_J(q", q"_, t; q^, q'_, 0)ps(?+, q'-, 0) , 



where 

J{q'l, q':, t- q'^, q'_, 0) = U(q':, q'^, t)W{q'l, q'_, t) (4) 

is the propagator with 

U(q':,q'^,t)^{q':m)\q'J (5) 

and 



ps{q+,q-) = {q+\ps\q-) 



(6) 



is the density matrix. 

According to Feynman lll4 ll5ll . the unitary propagator can 
be expressed as a path integral 



U(q'l,q'^,t) 



f 



Dq± exp 



TSs(q±,t) 

h 



(7) 



running over all paths which satisfy the boundary conditions 
q±(0) = ^±and^±(0 = ^". 



Ss[q,t] 



I dt'Lsiq,q,t') 
Jo 



(8) 



is the action associated with the path q{t') where Ls{q,q,t') is 
the Lagrangian describing degree of freedom of the system. 

As a system degree of freedom, we will specifically consider 
a harmonic oscillator of mass m and frequency loq. Then, the 
exact propagator takes the form 



, 1 ^^^s ^ 



exp 



jSf{q:,q:,t) 



(9) 



The classical action 



2 sin(woO 



(q'^^ + q'i) cos (woO - 2q'^q'^] 



(10) 

is completely determined by solutions of the classical equation 
of motion 

(11) 



qf = -mojlqf , 



while the fluctuations contribute only to the prefactor which is 
independent of q'^ and q'j_. 



For the harmonic oscillator in real space, we can thus con- 
clude that the propagator for the density matrix is determined 
by two independent solutions of the classical equation of mo- 
tion. The question now arises what can be said about the rele- 
vant trajectories in phase space where instead of two initial and 
two final points in position space one has only one initial and 
one final point in phase space, each comprising a position and 
a momentum. 

Among the infinitely many possible phase-space representa- 
tions il6ll we choose the Wigner function, which for our one- 
dimensional system degree of freedom is given by the Weyl- 
ordered transform of the density operator il6ll . 



Ws(r) 
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J 2nh 



expl-^M' U| 



PS 



q 

^"2 



(12) 



(3) Here, r - (p, q) denotes a vector in two-dimensional phase 



space. 

By means of the transformation (fT2l i one obtains from (O for 
the time evolution of the Wigner function 



Ws(r 



",0 = /dV 



Gw(r",r',OWs(r',0). 



(13) 



Introducing the difference coordinate q - q+ - q^ and the sum 
coordinate q - (q+ + q-)/2, the Wigner propagator appearing 
here as an integral kernel can be written in terms of a double 
Fourier transform of the propagator of the density matrix (HI 
along the difference coordinates 



Gw(r",r',f) = ^ I dq'dq"exp 



-{p'q' - p"q") 



, ,, q , q 

^U\q +—,q' + -,t 
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(14) 



A similar expression can be obtained with the position-space 
propagator replaced by the Weyl transform of the time evolution 
operator Onn . 



Gw(r' 



' 2nhJ 



d r exp 
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(15) 



where, analogous to ( fT2l l 



2.;.^'^Pl"W^^I 



U(t) 



q- 



(16) 



In this way, the classical trajectories q'^ determining the prop- 
agator (|9]l in position space become the phase-space trajecto- 
ries r± shown in Fig. [T] All such pairs share the same mean 
r = (r_ + r+)/2. For the harmonic oscillator, the phase-space 
difference vector f = r+ - r_ appears only linearly in the dy- 
namical phase. Thus, integration over f in ( fTsT i directly results 
in the exact Wigner propagator 



Gw(r",r',f) = 5[r"-rV,0]- 



(17) 




Figure 1 : For an undamped harmonic oscillator, the pairs of phase-space tra- 
jectories r+ as well as their sum trajectory r satisfy the classical equations of 
motion. r^,r' and r".r" denote the initial and final phase-space points of the 
trajectories, respectively, and f = r-,- - r_. 



can be exactly expressed in terms of the classical phase-space 
trajectories K - {r, R) of system and heat bath. They are ob- 
tained even in the continuum limit and for a thermal initial state 
of the bath by transforming to normal modes of the total sys- 
tem 12011 . It is only after tracing out the heat bath that quantum 
effects come into play in the time evolution of the Wigner func- 
tion. Instead of carrying out this analysis in phase space, we 
will rather follow the strategy presented in Sect. |2] for the un- 
damped case where we started from the position-space descrip- 
tion. 

For a dissipative system, one is usually not interested in the 
full dynamics but only in the reduced dynamics of the system 
degree of freedom. In order to integrate out the bath degrees 
of freedom, one needs to specify the initial density matrix of 
system and bath. For simplicity, in the main part of the paper we 
will restrict ourselves to factorizing initial conditions Ell 12211 . 
Then, the initial density matrix 



This result implies that each initial point r' in phase space is 
propagated along the solution r'^'(r', f) of the classica l eq uation 
of motion. All quantum effects, including tunneling llisll . must 
therefore be contained already in nonclassical features, such as 
an irreducible energy spread, of the initial Wigner function it- 
self. 

As we shall see in Sect. |5] the relevant trajectories will be 
different for the harmonic oscillator in the presence of damping. 
While the mean of a trajectory pair still satisfies the classical 
equation of motion, the two trajectories themselves will behave 
nonclassically and separate exponentially fast in time. How this 
comes about will be explained in the following sections. 

3. Influence functional theory 

We now turn the so far isolated quantum system S into a dis- 
sipative quantum system by coupling it to a heat bath character- 
ized by a temperature T . The complete Hamiltonian then is of 
the general form 



H - Hs + Hs + HsE 



(18) 



where the second and third term describe the bath Hamiltonian 
and the system-bath coupling, respectively. 

As long as the complete Hilbert space is retained, the evo- 
lution of the density matrix is still of the form (|2]i where the 
system density operator ps is replaced by the full density op- 
erator p. Correspondingly, in the time evolution operator ([T]i 
the system Hamiltonian has to be replaced by the full Hamil- 
tonian ( fTSl l. In position representation, ^ still holds with the 
appropriate replacements. In particular the coordinates are now 
<3 = (?> Q} and comprise the system coordinate q as well as the 
vector of bath coordinates Q. 

If one considers a harmonic oscillator bilinearly coupled to 
a bath of oscillators 11911 . a situation which we will discuss in 
detail below, the Wigner propagator 



p(0)=ps(0) 



exp(-;8//B) 



(20) 



is given by the product of the initial density matrix ps(0) of the 
system and the thermal density matrix of the heat bath at inverse 
temperature yS = l/k^T. Zb denotes the partition function of 
the bath. As we will show in Appendix A our reasoning can 



readily be generalized to a certain class of nonfactorizing initial 
conditions. 

The time evolution of the initial state ( l20b is obtained as a 
generalization of the considerations in Sect. [2] by substituting 
the single system degree of freedom by the ensemble of system 
and bath degrees of freedom. In order to obtain the reduced dy- 
namics of the system, one needs to trace out the environmental 
degrees of freedom. This can be done analytically if the heat 
bath is modelled by a set of harmonic oscillators with masses 
nij and frequencies aij whose coordinates are bilinearly coupled 
to the system coordinate. The Hamiltonian ( fTSl l then consists of 
the three contributions 123 1 



Hs = ^+ V(q) , 
Zm 



^B-Z 



;=i 



^ Pj + Irnjo^Q] 



2m 






(21) 
(22) 

(23) 



where Cj are the coupling constants. The second term in ( |23] ) 
corrects a potential renormalization induced by the coupling of 
the system to the heat bath. It turns out that the microscopic 
details of the heat bath and its coupling to the system appear in 
the reduced system dynamics only through the spectral density 
of bath oscillators 



/(w) = TT > 5(a) - Dj) 



j=i 



ImjOJj 



(24) 



Gv/CR", fi', 0-6 1'R" - 'R!'\'R', t)] . (19) For later convenience, we introduce two quantities which de- 



pend on this spectral density. The friction kernel is defined as 

77(f) = - I dw— - cos(wf) (25) 

TT Jo w 



and the correlation function of the noise 112411 induced by the 
coupling to the heat bath is given by 



K{z) 



dw cosh[w(ifi,8-iz)] 



— I{0J) 

Jo ^ 



sinh(fij6w/2) 



(26) 



where z is a generally complex time. 

Tracing out the heat bath, one finds for the time evolution 
of the reduced density matrix of the system starting from the 
factorizing initial state ( l20l l 112 U 12211 



P&{q'i, q'-, t) = I Aq'^ dq'_J(q", q", t; q'^, q'_)ps(q'+, q'^, 0) . 

(27) 
Here we have introduced the propagating function which can 
be expressed in terms of a path integral over the system degree 
of freedom as 



J(q':,q';,t;q'^,q'_) ^ ^ j Dq^Dq^ 



X exp ^ - [S s(?+) -Ssiq-)]} Tiq^, q-) . 



(28) 



The partition function Z appearing here is an effective partition 
function of the damped system defined as the ratio of the par- 
tition functions of system plus bath and of the heat bath alone. 
The influence of the heat bath in ( l28l l is contained in the influ- 
ence functional 



T{q+,q-) = exp 
with the exponent 

^[q+ 



--0(q+,q-) 



(29) 



q-] - I ds I duK'{s -u)[q-^(s) -q^(s)] 
Jo Jo 

X [q+(u) - q-{u)] 

+ - I ds I duj]{s-u)[q+(s)-q4s)] (30) 
^ Jo Jo 

X [q+(u) + q-(u)] 
+ ^[q'++q-]j ds[q^{s)-q-is)]. 



The prime in the first line indicates that the real part of the noise 
correlation function (|26]| should be taken. 



4. Quantum damped harmonic oscillator 

While the results reviewed in the previous section are valid 
for a general system degree of freedom, we will now specifi- 
cally consider a damped harmonic oscillator with 



2 "2 



V{q) = -(^Qq 



(31) 



In this way, we will be able to generalize the considerations of 
Sect. |2] concerning the relevant trajectories in phase space. At 
the same time, this prevents making direct use of semiclassical 
phase-space dynamics based on the van-Vleck approximation 
as in Refs. flldlllll . since this only applies to sufficiently anhar- 
monic potentials. 

As for the propagator in the undamped case, the path-integral 
expression for the propagating function ( l28l l is evaluated by 
an expansion around the paths maximizing the complex action. 
The dependence on the initial and final coordinates is entirely 
determined by these paths while the fluctuations only yield a 
time-dependent prefactor. For a harmonic oscillator, the com- 
plex action in ( l28l l is stationary for trajectories satisfying 

1 d r' 

mqj,(s) + mcolqj,{s) + -— I duj](s - u)[q+(u) - q^iu)] 

1 d r 

+ 23~l duT]{s-u)[q+{u) + q-{u)\ = 



i I duK'{s — u)[q+{u) — q-{u)] . 
Jo 



(32) 

As in the undamped case, the paths are subject to the boundary 
conditions ^'±(0) = q'j_ and ^±(0 - q'l- 

The two equations of motion (|32] | replace the equations of 
motion ( fTTT i obtained for the undamped case. As a consequence 
of the coupling to the heat bath, the equations of motion ( l32l l 
now include damping terms depending on the friction kernel 
( |25] |. In addition, an imaginary nonlocal force appears indicat- 
ing the occurrence of decoherence. For linear systems, it turns 
out that the imaginary part of the trajectories does not need to 
be considered and that their real part is sufficient to obtain the 
propagating function 125112611 . As we shall see in the sequel, 
neither of the two paths follows a classical equation of motion 
and their separation grows exponentially fast. This somewhat 
surprising behavior is a consequence of the coupling to the heat 
bath. 

In order to make the discussion transparent, we now consider 
the special case of Ohmic damping in addition to the assump- 
tion of factorizing initial conditions and the restriction to the 
real part of the equations of motion. We thus set /(w) - myo) 
and y{t) -rj(t)lm - 2y6(t), so that the equations of motion (l32l l 
reduce to ll22ll 

q± + Ci>lq± +yq+^0, (33) 

where the damping term couples the trajectories q+ and q^. It is 
interesting to note that yq+ acts actually as a driving instead of 
a damping in the sense that the separation between trajectories 
grows exponentially (see Fig.|2]i. This can be seen more clearly 
by decoupling the two equations of motion, using sum, q - 
{q+ + q-)/2, and diff'erence coordinates, q = q+ - q^. The 
equations (l33T l then read 



q + yq + cij^q = 
q - yq + ojf^q = . 



(34) 



As in the undamped case, the sum trajectory corresponding 
to the paths q+ and q^ obeys the classical equation of motion. 



which here takes a time-local form because we have assumed 
Ohmic damping. In contrast to the sum coordinate q which de- 
creases exponentially in time, the difference coordinate q grows 
exponentially so that we obtain a hyperbolic dynamics in the 
{q, §)-plane. As a consequence, the trajectories q± do not obey 
the classical damped equation of motion. The solutions of the 
equations of motion ( l34l i read 



, GAt-s) , „G+W 



G-(t) 



GAt) ' 






(35) 



G+(0 



G-(0 ' 



where 



G±{f) = — exp(+^f) sin(wdO (36) 

Wd \ 2 / 



and 6J^ = Wq - y^/4. We remark that by choosing appropri- 
ate functions G{f), more general linear damped system like the 
parametrically driven damped harmonic oscillator Ii27i - t29ll can 
be studied with solutions of the form (l35T l. 

We now return to the propagating function which was given 
in ( l28T l in terms of a twofold path integral. It is instructive to 
decompose the exponent into two parts 



S(q",q",t;q',q') - S [ + S2, 



where 








Si = m 


(,'?., ■■f.if 


-"" GM) 


'''C.(0 



(37) 



(38) 



is obtained by evaluating the action of the system degree of 
freedom along the trajectories given by (l35T l while 



'.i:^i 



du K'{s — u)q{s)q(u) 



(39) 



arises from the influence functional ( 1291 ), i.e. by the interac- 
tion of paths at different times through the coupling to the en- 
vironment. The significance of this decomposition will become 
clear in the following section where we discuss the results of 
the present section from a phase-space point of view. 

5. Quantum damped harmonic oscillator in phase space 

The undamped time evolution ( fT3b of the Wigner function 
can immediately be transferred to the dissipative case if we re- 
late the Wigner propagating function to the propagating func- 
tion by means of 



Gw(r",r',0 



2^ J'^^' 



dq" exp 



-Xp'q -p"q") 



(40) 



xJ(q",q",t;q,q'), 



where we recall that r = {p, q) is the phase-space vector. In- 
serting Q valid for the unitary case into (|40] | one recovers ex- 
pression (fT4l l for the Wigner propagator We continue to re- 
strict ourselves to factorizing initial conditions but present in 



Appendix A the generalization of ( |40] | to nonfactorizing initial 
conditions. 

Before analyzing the Wigner propagating function for Ohmic 
damping, we discuss the decomposition ( l37T i of the exponent of 
the propagating function. The contribution ( l38l l is linear in the 
difference coordinates q' and q" . Performing the transforma- 
tion ( l40l i, we therefore arrive at the Wigner propagating func- 
tion 

Gw(r", r', f)^6 [r" - r'='(r', f)] , (41) 

where the classical phase-space trajectory 



p%) = G^{t)p' + m 



1 



{GAt)f 

G+(0 G_(0 



q%)^^^p'+G,{t)q', 



(42) 



with G±{t) defined by (|36] | is now damped. While (l42l i satis- 
fies ^''^'(0) = q' as expected, the initial momentum is given by 
p'^'(O) - p' - myq' . This initial slip is typical for factorizing 
initial conditions ll26ll30ll . 

Employing the Wigner propagating function (1411 1 amounts to 
adding a velocity-dependent force in the system Hamiltonian 
as was proposed by Caldirola and Kanai (see e.g. the review 
1 3l|] for an account of this kind of phenomenological approach) . 
However, the Wigner propagating function (l4TT l accounts only 
for part of the exponent of the propagating function. The second 
contribution ( [39] l is quadratic in the difference coordinate and 
limits their contributions. As a result, the delta function in (|4TI) 
will be broadened into a Gaussian 



Gw(r",f;r') 



^^'^p{-2^['-"-'''«r^['-"-'-"^^^]}' ^"'^ 

whose center moves along the damped classical trajectory (l42T l. 
The matrix appearing in the exponent is given by its compo- 
nents 

111 =a(f) 
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GAt) 

G+(0 


2nhK(tyi^ 



Z12 = I21 = -m§4? [a{t) + bit)] 



I22 - m 



2[G+(0]' 



GAt) 
[a{t) + 2b{t) + c{t)] 



(44) 



[GAt)]'' 
and A(f) = det(Z)/m^ - a{t)c{t) - b{tf. The functions 

a{t)^[G^{t)f^{t,t) 



b(t) = G+(OG+(f) 



d'Vit, t') 



dt 



I' A 



(45) 



c{t) = [G+(0] 



2 d'^'Vit, f) 



dtdf 



f/t 



can all be expressed in terms of a single function 

G+(t - s) G+{t' - u) 



"Vit, t 
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Au K'{s — u)- 



G4t) GAt') 



(46) 



This function is completely determined by the thermal position 
autocorrelation function {q{t)q{0)) and its time derivatives. Ex- 



plicit expressions are given in Appendix B and the interested 



reader may also want to consult Ref. 112611 for further details. 

The Gaussian form of the Wigner propagating function ( |43] ) 
is a consequence of the linearity of the harmonic oscillator 
damped by the coupling to a harmonic heat bath. A similar ex- 
pression has therefore been found in the Markovian limit ll32ll . 



Similarly, as indicated in Appendix A the result can be gen- 
eralized to the case of nonfactorizing initial conditions. Again 
one would find a Gaussian Wigner propagating function. 

From ( l40l i it follows that pairs of trajectories q-t satisfying 
the equations of motion (l33l l and leading to a sum trajectory 
connecting the initial phase-space point r' to the end point r" 
contribute with a weight determined by ( [39] ). As discussed in 
the previous section, according to ( l33l l these pairs do not obey 
the classical damped equation of motion. The same is true in 
phase space. The hyperbolic character of the solutions of ( l34b 
already contains all quantum effects that arise upon tracing out 
the heat bath. The equations of motion ( l33T l can therefore be 
lifted to the phase space of the central oscillator by defining 
P± = q±/m, 

p^ = -ma)lq± - myq^ 
. P. (47) 

m 
so that the sum trajectory follows the classical equation of mo- 
tion of the damped harmonic oscillator This result was also 
found by Ozorio and Brodier lll2ll for the Markovian case de- 
rived directly from the Lindblad master equation. 

In Fig. |2]we depict the time evolution in phase space of two 
trajectories q± indicated by + and - together with the corre- 
sponding sum trajectory shown in black. Due to the damping, 
the sum trajectory for long times approaches the origin of phase 
space. The trajectories q± grow exponentially for long times 
and therefore clearly behave nonclassicaUy. Although in Fig. |2] 
the paths q+ and q^ have started on the same side of the origin 
of phase space, for long times they are found opposite to each 
other This is a consequence of their exponential growth and 
of the fact that their midpoint approaches the origin of phase 
space. 

We close our discussion of the phase-space properties of the 
damped harmonic oscillator by considering how the thermal 
equilibrium state is approached for long times. First, we notice 
that in the Wigner propagator (l43l l the dependence on the initial 
phase-space coordinates r' disappears in that limit because the 
sum coordinate then approaches the origin of phase space. Fur- 
thermore, the long-time behavior of the functions (l45T l is given 
by 



a{t) 



bit) 
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where {q^} and {p^} are the second moments of position and 
momentum, respectively, in thermal equilibrium. Inserting 




Figure 2: The time evolution of a pair of phase-space trajectories r+ marked by 
+ (depicted in blue) and by — (depicted in red) is shown together with the cor- 
responding classical sum trajectory depicted in black. While the sum trajectory 
decays to zero for long times, the trajectoiies r+ grow exponentially. 




Figure 3: Isosurface of the time-dependent Wigner propagating function (43) 
for y/ojo = 0.3 and k^T = StiojQ. Position and momentum are scaled with 
respect to the square root of the respective equihbrium second moments. 



these expressions into ( |43] ), one obtains the thermal Wigner 
function of the damped harmonic oscillator 



Wi,{p,q)^ 



1 



■ exp 



q 



(49) 



2;r«^2></,2))i/2 ^\ 2{p^} 2{q^} ) 

In Fig. |3]we illustrate the time evolution of the Wigner prop- 
agating function for y - 0.3(^0 and k-^T - Sftwo by means of an 
isosurface. The function ( |46] ) has been evaluated with the high- 
temperature approximation K'{t) - {2y/hp)S{t). The propagat- 
ing function evolves from an initial delta function to the thermal 
Wigner function ( |49] ) which has a circular cross section because 
position and momentum are scaled with the square roots of the 
respective equilibrium second moments. 

6. Conclusion 

Our analysis of the propagating function of a damped har- 
monic oscillator in phase space has shown that pairs of trajec- 



tories contribute which, in contrast to the undamped case, are 
nonclassical in nature. In the course of time the distance be- 
tween the two trajectories of a pair grows exponentially. This 
behavior is typical also for nonlinear systems 11311 so that the 
damped harmonic oscillator can serve as an exactly solvable 
reference case. The harmonic oscillator coupled to a linear heat 
bath is however special insofar as the linearity of the system im- 
plies that the sum of the pairs of trajectories satisfies the classi- 
cal damped equation of motion. 

Furthermore, it was shown that it is the nonlocal interac- 
tion ( l39l l appearing in a reduced description of the system de- 
gree of freedom which is responsible for the broadening of 
the delta-like Wigner propagator found for the undamped case 
into the Gaussian propagating function. While most of the ex- 
plicit results have been presented for the special case of fac- 
torizing initial conditions, the generalization to a broader class 
of initial states is straightforward along the lines presented in 



Appendix A and the Gaussian nature of the propagating func- 
tion remains unchanged. 
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Appendix A. Propagating function for nonfactorizing ini- 
tial conditions 



In the main part of this paper, we have restricted our con- 
siderations to the case of factorizing initial conditions for the 
sake of simplicity. The generalization to nonfactorizing initial 
conditions is a bit more tedious but straightforward. In this ap- 
pendix, we briefly derive the phase-space representation of the 
time evolution of a nonfactorizing initial Wigner function Ii33i1 . 

Specifically, we will consider nonfactorizing initial condi- 
tions where system and bath together are in their equilibrium 
state thus accounting for initial correlations. Then, operators 
taken from the Hilbert space of the system are allowed to act 
in order to generate a nonthermal initial state for the system. 
On the formal side, this preparation in comparison with factor- 
izing initial conditions has two consequences. Firstly, in ad- 
dition to the two real-time paths q±, an imaginary-time path q 
appears. Secondly, a preparation function A{q',q',q,q) joins 
the real-time and imaginary-time paths. Here, q and q refer to 
sum and difference coordinates of the imaginary time path. The 
preparation function function describes the initial preparation 
and depends on the matrix elements of the system operators. 
For more details, we refer the reader to Ref. i26ll. 



The generalization of (l27T i to the nonfactorizing initial con- 
ditions just described reads 



ps(q'\q'\t) 



= fd^'d^' 



dqdq J(q", q" , t\ q ,q' , q, q) 



(A.l) 



X A(q ,q ,q,q) 



Now, the propagating function also depends on the endpoints of 
the imaginary-time path. Introducing the Wigner transform of 
the preparation function 



A^ 



^{p',q',p,q) = —-—2 I d^'d^exp 

X A{q' ,q' ,q,q) , 



-ipq-p'q') 



(A.2) 



we obtain for the time evolution of the Wigner function after 
carrying out the Fourier transform with respect to q" 



Ws(p",q")= I dp'dq'dpdqGw(p",q",t;p',q',p, 



q) 



(A.3) 



xAw(p',q',p,q). 



By comparison with ( lA.ll ) one finds for the relation between the 
propagating function introduced in ( lA.ll ) and its Wigner trans- 
form 

Gw(p",q"j;p\q\p,q)^ —- I dq"dq'dq 

i , .^. 
X exp 



-(pq -p q -pq) 



J{q ,q ,t;q ,q ,q,q) . 

(A.4) 
For the special case of factorizing initial conditions, the coordi- 
nates q, q and the momentum q are to be disregarded and one 
arrives at the relation ( l40b between the propagating functions in 
position and phase space. 



Appendix B. Propagating function and position autocorre- 
lation function 

The Gaussian nature of a harmonic oscillator coupled lin- 
early to a bath of harmonic oscillators implies that its reduced 
dynamics can be expressed completely in terms of the thermal 
position autocorrelation function 13413511 



C{t) = {q(t)qm 



-r 



dot 



S{t) + iA(t) 

cjy(—icbi) 



exp(-ia>t) 



(co^ 



■ Wq)2 + CLi-y{-ia>)^ 1 - e\p{-/3ha)) ' 

(B.l) 
S (t) and A{t) denote the symmetrized and antisymmetrized cor- 
relation functions and correspond to the real and imaginary part 
of C{t), respectively. y(z) is the Laplace transform of the fric- 
tion kernel dZSl l divided by the oscillator mass m. The anti- 
symmetric correlation function is related to the function G+{t) 
introduced for the special case of Ohmic damping in ( |36] ) by 






(B.2) 



where 0(f) is the unit step function. The second moments of 
position and momentum appearing in ( |49l ) are related to the 
symmetrized correlation function by {q^) - 5(0) and {p^) = 
-mS{0), respectively. For the latter to be finite, the Laplace 
transform y{z) requires a high-frequency cutoff. 
The functions ( l45T l are found to read 12611 






^_lS{t)f 



h G^t) 






+ ^[G+(0]' + 25(0G+(0 



{q^y~ 



!§.<;,„)*„) 



-«,«.«. ^} 



«'>-ii^"'' 



'^« = ir|<^>[G;(oF"^-<?) 

(B.3) 
If one takes into account that G+{t), S{t), and their derivatives 
decay to zero for long times but not the ratio G+{t)/G+{t), one 
immediately finds the asymptotic expressions (l48b . 
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